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Abstract 

A magnetic monopole is placed at the centre of a ball whose surface 
^ ' S^ is tiled by the symmetry group, F, of a regular solid. The quantum 

\Q I mechanics on the two-dimensional quotient S^/F is developed and the 



O 



monopole charge is found to be quantised in an expected manner. The 



O ' heat-kernel and C^functions are evaluated and the Casimir energy is 

Q\ , computed. Numerical approaches to the calculation of the derivative 

of the Barnes C^function are presented. 
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1. Introduction 

In previous work, [1] , we tiave discussed quantum field theory on orbifold factors 
of spheres S"/r where F is a discrete subgroup of the orthogonal group, 0(n + 1). 
In this paper we wish to present an extension in the two-dimensional case (n = 2) 
to the situation where there is a uniform radial magnetic field passing through 
the surface of the sphere. This can be thought of as due to a magnetic monopole 
at the centre of a ball in an embedding R^. The motivation is partly to investi- 
gate the interplay between a magnetic field and a non-trivial topology/geometry 
induced by identification. Specifically we would be interested in what happens to 
the Dirac quantization in topologically interesting or singular manifolds (orbifolds). 
There is also continuing statistical mechanical interest in magnetic fields and two- 
dimensional domains. 

The quotient group is the complete symmetry group of a regular solid and can 
be generated by refiections in the three (concurrent) planes of symmetry of the solid. 
The elements fall into two sets depending on whether they contain an even or an 
odd number of relections. The even subset forms the rotational subgroup denoted 
now by r G SO (3). The odd subset is denoted by Fp If 7' is any fixed element 
of F' then as 7 runs over F, 77' exhausts F'^^. In particular we can choose 7' = cr 
where cr is a refiection in a symmetry plane and so the complete group is 

F' = F U F; = F U Fa. (1) 

The standard classification of finite subgroups (reflection groups) is given by 
Meyer, [2], and we use his notation. The general construction of the quotients, 
S^/F', which in this case are certain geodesic triangles on S^ is sketched later. 
The vertices of these triangles are singular points. Joining them to the origin of 
the ambient M produces singular strings and our analysis can be extended to this 
three-dimensional setting [3]. The two-dimensional models that we study can be 
thought of as toy models for more general 'textures' in higher dimensions as laid 
down by Kibble, [see e.g. [4]). 

Our main interest is in setting up the general framework and then applying it 
to some specific calculations such as the evaluation of vacuum energies. This will 
involve a certain amount of technical manipulation, especially with the properties 
of the Barnes (^-function. This function appears quite commonly, particularly in 
spherically symmetric situations and in the presence of magnetic fields or harmonic 
oscillators and we expect that our methods will have an applicability beyond the 
immediate one here. 
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A very brief summary of our findings is given in section 12. 



2. Modes and group actions on the full sphere. 

As modes on the full sphere, we can take the angular part of the Schrodinger 
equation solutions derived long ago by Tamm [5] and Fierz [6] . It would be more rig- 
orous to use a fibre bundle formulation (Greub and Petry [7] , Wu and Yang [8]) or 
geometric quantisation [9,10] but this is unnecessary for our purposes. For definite- 
ness, we employ the modes denoted by (Y^^)^ in Wu and Yang, corresponding to the 
string running down the negative z-axis. The modes are, up to normalisation, the 
SU(2) representation matrices Vl^*_q{(p, 9, —(p) with / = |g|, |g + 1| . . . , —l<m<l 
[11]). 2q is the monopole number with q = eg/h, and 2q G Z. If the string runs up 
the positive z-axis, the modes are r'^\lq((/), 6*, (p). 

The corresponding eigenvalues of Hg2, the angular part of —(V — zeA)^, are 

Xi = l{l + l)-q'=(l + l^ -\-q' (2) 

with degeneracy 2/ + 1. 

It is helpful to give the explicit form of the angular eigenfunctions in spherical 
polar coordinates, 

n^^(^,</')=A^,.^sinl^+'-l(0/2)cosl^— l(0/2)pl!+7|i^-7^^^^ 

(3) 
where N i^ is a normalisation constant and the P^'^{x) are Jacobi polynomials. In 
the Wu-Yang formalism these are the solutions (sections) in the upper hemisphere. 
In the lower hemisphere the potential is taken to have a string along the positive 
z-axis. The two sets of solutions are related on the equator by the factor exp{i2q(j)). 
Making this single valued gives the quantisation condition on q, in this approach. 

The basic means of finding the modes is to write -ffg^V' = -^V' explicitly as a 
differential equation in spherical polar coordinates, 

1 d . „dih 1 d'^th . 2(7 dth 9I — cos6' , 

sm^-^ + — ^-— J+z- ^-—^-q^- ^^ = A^, (4) 



sin^a^ de sin^edcj)'^ l + cosed(j) ^ 1 + cos^ 

and solve it assuming the separation 

1^(9, (/)) = exp(iu(j))P{cos9). Using this method we can suppose initially that q 
is arbitrary which leads to the same solutions as in (3) but which are characterised 



by two integers u G Z"*" and f G Z, [5]. The relationship to the SU{2) labels I and 
m is given by 

I = u+ -[\q + m\ + \q — m\), v = m — q (5) 

and thus it would appear that we could have any value of q as long as u and v are 
integers (with suitable I, m &M.). However, in this case the solutions would not be 
single-valued and all the solutions would vanish on the string axis. Setting 2(j G Z 
gives the same values of / and m as before. 

As usual, when one tries to adapt wave-functions to some symmetry (here V) 
there is the problem that the magnetic potential. A, might not possess the same 
symmetry so that compensating gauge transformations are necessary. Peierls, [12], 
calls this process 'umeichen'. It is a well known situation, with extensive discus- 
sion, [13], which we have encountered in an earlier calculation on the tetrahedron, 
[14], an orbifold factor of the plane. A similar procedure has also been applied to 
factors of the Poincare half-plane. 

In the present case, under the action of F', the string is rotated and reflected 
and has to be brought back to its original position if we are to implement the iden- 
tification in S^/r' consistently. The equation that contains the relevant information 
is the behaviour of the modes under arbitrary rotations-reflections. For the moment 
we deal with the easier, pure rotational case, g G S0(3). Then the behaviour is an 
elementary consequence of the SU(2) group combination law and one has explicitly 
(Wu and Yang [8], Frenkel and Hrasko [15]) 

YiZi9-'r) = e^^^(^) y«(r) vl^^,{g), g G S0(3) . (6) 

The exponential factor is the compensating gauge rotation with 

Ag{?) = a{g)+^ig)-ng (7) 

where a, (3 and 7 are the Euler angles of the rotation g and O^ is the solid angle 
subtended by the geodesic triangle on the (unit) sphere cut out by — r, the string 
axis, n, (here the negative ^-axis) and the rotated string axis, gn, 

^s(r) = ^(fi-n, n;r), r = (r, r). 

(Our conventions regarding rotations and actions are generally those of Brink and 
Satchler [16].)The gradient of O gives the gauge transformation between the poten- 
tials of the two strings. 



An alternative expression for A is (Wu and Yang [8]) 

Kg{v)=(t)g-(t)-Ag (8) 

where A is the angle at — r in the above mentioned triangle on S^ and g~^r = 

In order that (6) be consistently iterated, it is necessary that the group com- 
bination (cocycle) condition, 

A^J?)=A^(?)+AJ<7-i?), (9) 

should hold, and this can be checked directly from (7). The geometrical details are 
to be found in the useful article by Frenkel and Hrasko [15]. 

The magnetic rotation operator, T , is defined, on scalars, by the action 

(T^^)(r) = e-^^^(^)V'(^-'r), (10) 

or, equivalent ly, 

(r|T^ = e-^^^(^)(<7-^r| . (11) 

In the present spherical case, because of (9), the magnetic rotations provide a true, 
as opposed to a ray, representation of the double of S0(3) (denoted SO* (3)) in the 
sense that 

and 

TE = {-ifn, (12) 

where E' is a 27r rotation. This is not unexpected considering that for q a half 
odd-integer all the monopole harmonics have half odd-integer angular momentum. 
SO* (3) is isomorphic to SU(2). 

Magnetic rotations on the plane [17] also provide true representations but not 
so translations, unless the fiux through a fundamental domain is quantised [18]. 

We record the action of T on the modes which is easily obtained from (10) 
and (6), 

r,^J^.(r) = y£(?)i^!:L,(^). (13) 



3. Parity and reflections. 

The singularity (string) preserving magnetic parity operator on monopole wave- 
functions is defined by the action, [15], 

(nV')^ = e-*^"^-"'"''"^ V'(-r) (14), 

under the inversion l : r ^ —r 

On the modes we have, as an easy calculation shows, 

nr« = (-i)'y«^. (15) 

The parity operator can be used to extend the magnetic rotation operator T 
to include reflections. The reflection, a, in the plane with normal t can be written 
as a rotation through tt about the axis t combined with the parity inversion l. This 
can be written a,s a = g^ i = cg^. The extension of T to reflections is thus defined 
by 



From (10) and (14) 



and on the modes. 



(T»(r) = e-^'^"^«+^'^- (R ^)(-r), (16) 



(^.C')(^) = (-1)' ^-U(^) ^w(^J- (17) 



4. Construction of the domains S^/r' 

It is convenient now to formalize briefly what we mean by the space S^/F' 
where F' is a finite subgroup of 0(3). A region of the sphere JF e S^ is called a 
fundamental domain for F if it satisfies the following criteria, 

(i) JF is open in S^ , 

(ii)^n7^ = 0_V7eF'-H}, 

(iii) S2 = U^er' 7^- 
We also assume that JF is connected. 

Physically, JF |J dJ-' represents the space on which our theory is defined, and all 
the copies of JF on S^ must be physically equivalent. 



The space S^/F' is the the sphere S^ with the points r and 'jr identified for all 
7 G r'. If r' acts freely on S^, i.e. there are no fixed points, then S^/F' is closed 
and can be taken to be JF with identified boundary. If there are fixed points, these 
will be contained in dJ-' and can be considered as boundary singular points of S^/F'. 

When F is purely rotational the fixed points form a discrete set of which there 
are two, or three, in dJ-'. For the extended, refiection groups, as already indicated, 
the boundary dJ-' is constructed from the intersections of two or three refiection 
planes with the sphere. Thus there are transformations that map JF into adjacent 
domains and which leave part of the boundary of JF fixed. The conclusion here 
is that, unlike the rotational case, the boundary dJ-' is a real boundary and the 
physical manifold is termed a Mobius triangle. 

5. Projection to S^/F 

Still restricting to purely rotational F, we will now formally project everything 
down to S^/F in a more or less standard fashion. 

In order that a function V'(r) on S^ project down consistently to S^/F, it is 
necessary that it satisfy the magnetic periodicity condition 

Tg.^{r)=a{g')^{r), geV, (18) 

where a{g*) forms a one-dimensional representation of F*, the double of F, and 
labels the projection. For ease we have put r in place of r. We show later that there 
is a minimal choice for a{g*). 

If q is integral, there is no need to introduce the group doubles. However, if 
one does, there is a simple duplication, which is easily dealt with in practice as we 
can show with the following, somewhat superfiuous, constructions. 

Since the action of F* covers the sphere S^ twice, we introduce the trivial double 
covering, S^ , of two identical copies of S^, with S^ /F* = S^/F. 

The modes themselves do not change sign on a 2n rotation, even when / is half 
odd-integral, and neither does the wavefunction. So i(^*{r) = t(j*{Er), where E is a, 
2n rotation, and therefore, il^{r) = V'*(r)- 

We then have, somewhat non-rigorously, for the function, ^(r), on S^ /F* the 
numerical equality 

-i/;(r) = V*(r) (19) 

where we do not distinguish between the projected and original coordinates, both 
being denoted by r. The r on the right belongs to the fundamental domain of F* 
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on S^ , which is, of course, isomorphic to that of F on S^. Then 

T^.V(r)=a((7*)^(r). (20) 

The heat- kernel on S^, written in mode form, 

oo 

i^S^(r,r';r) = J^ e-^'^TY [Y«(r)Y«t(r')] , (21) 

propagates arbitrary wavefunctions on S^ and satisfies the important switching 
relation 

Ks (r, g-h'; r) = e'^'^^.^'-'^-^'^^- « K^ {gr, r'; r) 

which realises in coordinate representation the operator rotational symmetry 

(see (11)). 

The heat-kernel that propagates wavefunctions on S^ /T* obeying (20) is, by 

general theory, 

Ks-/r.(r,r') = J] a{g-)T^,-^ Ks..{r,r') (22) 

ger' 

which is referred to as the pre-image form of this propagator in terms of that on 

The double group F* can be decomposed 

F* - F U ^F 
and the sum over F* reduced to a sum over the subset F 

i^s-/r. (r, r') = Y, Hd) T,-i + a{Eg) T^Eg)~^) Kg- (r, r'). (23) 

We know, (12), that T^ = {-ifn and so from (20), a{E) = {-if'^. Hence 
a{E) T^ — 1 and so 

i^s-/r.(r,r') =2^a(<7)T^-ii^s-(r,r'), (24) 

ger 

which means we can finally write what was almost obvious from the start, the 

preimage sum, 

i^svr(r, r') = J^ a{g) T^-. K^. (r, r'). (25) 

gev 
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The factor of 2 is a normalization (volume) factor between S^ and S^, i^s2» = '2.Kq2. 

This result shows that we can effectively ignore the complication of the double 
group and just proceed with V as usual. 

Given an arbitrary field (/)(r) on S^, a quasi-periodic field on S^ and hence, via 
(19), a field on S^/F is constructed by the projection 

i^{^) = ^.Y.''^3)Tg-.4>{^). (26) 

In particular the projected {i.e. adapted) modes are 

y;m(r) = -^Y.a(g) (r,-,ym)(r) = VV\Y^^\^)P^'^ (27) 



rl.er 



where, using (13), 



p 

rn'm ipi 



ger 

is a hermitian projection matrix, P^ = P. Since the eigenvalues of P are and 1, 
this shows that in a suitable, i.e. diagonal, basis the projected eigenf unctions are 
a subset of the unprojected ones, a general result and independent of the magnetic 
field. 

Depending on whether q and / are integral or half odd-integral together, the 
representations a and V of the doubled elements, £^7, have the same or opposite 
sign as those of 7 and we can see explicitly that it is adequate to restrict the sum 
in (28) to r, as we have done. 

Let us pursue this a little further to make sure everything is satisfactory. It 
is much more elegant to express everything in abstract form but we will retain the 
coordinate representation. The discussion is a textbook one in applied group theory. 

From (27, using periodicity, one easily derives the integral 






V 

Diagonalising A, A = UDU~^^ introduces the linear combinations Y^ '>U^ 
which, from (27) equals \/\r\YUD showing that certain linear combinations of 
the modes on S^ vanish and others do not, after adaptation which takes on the 
nature of a filtering process. We write for the nonzero combinations 

^ qa V I I qm^ ma, \ ^1 



which form a complete, orthogonal set on S^/F. 

The factors of ^/\^\ have been chosen to give suitably normalised eigenfunctions 
on SVr so that 

i^svr(r,r';r) = J^ ^""'^^ [3^i'Hr)3^«t(r')] . (30) 

i=\<i\ 

The matrix U reduces the /-representation of F* into irreducible ones and, from 
(28), the range of a, i.e. the degeneracy, is the number of times the irreducible a- 
representation occurs in this decomposition. 

5. Effect of fixed points 

The standard theory of coverings applies when the covering group acts freely. 
This is not so in the present case and the result is a singular space - an orbifold. 
We can still maintain the standard covering terminology, as we already have done 
when referring to fundamental domains, in an obvious way by firstly removing the 
fixed points, so that the standard theory applies, and then extending the action of 
the covering group, F, to these points. 

We now investigate the implications of the periodicity condition (18) when 
extended to the fixed point set of F. For any 7 G F there are two fixed points on the 
sphere, ±r . Going back to the mathematical result (6), we see that the monopole 
modes on the full sphere satisfy the identity 

y,ti±r,) = e^^^^^^'-^) l^£(±r,) I^^L(7), (31) 

where 7 is the rotation about ±r^, obviously. 

Using the relation (29) we can determine the action of the magnetic rotation 
operator T^ on the 3^, 

T^y{r) = ^/\r\{T^Y){r)U = e-"^^-'^''^ ^/\^\Y{-f-^r)U = e-"i^^^''^y{-f-^r) 

as expected. However we also have, from (27), 

T^yir) = a(7)3^(r) 

and therefore 

e-^^^^('-)3^(7-V) = a(7)3^(r) 
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and of course 

If there were no fixed points there would be no problem to discuss. However if 
r = r^ is fixed by 7 we have the consistency condition 

This says that the phase is undetermined at the fixed points. If the fixed point 
is removed there is no problem. We are then effectively working arbitrarily close 
to the fixed point and all the points 7r, 7 G F, belong to different images of the 
fundamental domain. When the fixed point is put back and the limit r — i> r taken, 
all these images coalesce and the phase becomes indeterminate, unless something 
special happens such as e~^'^ ^('■7-* equalling 0(7) or if the wave function has a node 
at r^. We analyse this situation further. 

Since F can be generated by a set of cyclic rotations, it is helpful firstly to take 
the case when F is cyclic, Cj^, about the ^-axis. Everything in (28) is then explicit. 
Let the generator of Cj^ be 7, (7^^ = 1). Also let the generator of the double group 
Z' be 7* with (7*)^ = E and {'j')'^'' = 1. E is the 27r doubling rotation. 

The representations a{g*) are 

a((7')P)=e2-W-, p = 0,l,...,2z/-l, (33) 

where, for integral q, the label r, is in the range to z/ — 1, while for half-odd integral 
q, r= (2s + l)/2 with < s < v - 1. 
^mm'(5'') is diagonal 

^^L ((7*)^) = e-^-^-/^^ 5_„ -/ <m,m'< /, (34) 

and -P^^/ is easily found, 

1 1 _ p2TTi{m+r) 

p = i___ A CiK) 

^ mm' 1/ I — ^2'Ki{m+r)/u ^mm'- V^^J 

This expression vanishes unless m + r is mod z/, in which case it equals unity 
making the filtering obvious. 

It is comforting to check things by looking at the explicit forms of some modes. 
For Q = 0, the modes Ym are ordinary spherical harmonics. At the north pole only 
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the 171 = component survives which then forces r to be zero and (32) is trivially 
satisfied since 0(7) = 1. 

One can proceed generally but let us use the modes exhibited in Wu and Yang 
section 7. For q = 1/2 we see that only the m = —1/2 survives at the north pole 
making r = 1/2 and now to check (32) we need the expression for A. For a rotation 
through uj about the ^-axis, from (7) or (8), A = — a;, and again the check works. 
For q = 1, the nonzero mode corresponds to r = 1. 

In general one finds that r = q and looking at the (p dependence of the modes, 
Yrhq, (3), i.e. exp (i(m+q')(/)) , we see that the modes Y^^'-^ are periodic as (p increases 
by 27t/u and so, therefore, is the wave function. 

It should be noted that if (8) is used, the azimuthal angle of the north pole 
changes by uj even though it is a fixed point. The string's location may be un- 
changed, but a nonzero compensating gauge transformation is still required. 

It is clear geometrically that the same results will hold if Cj^ is cyclic about 
any axis so that it is consistent to set 

a(7)=e-*'?^^('"^\ (36) 

for any F. The consequence of replacing r by — r is discussed in the next section. 
The minimal choice (36) clearly corresponds to untwisted fields in the sense 
that the more general form 

a(7) = e-^^^^^'"^U(7), (37) 

where 6(7) is some nontrivial representation of F*, encodes physics over and above 
the magnetic monopole, such as Aharonov-Bohm fiuxes through the fixed points. 
In this latter case, (18) implies that the wavefunction vanishes at the fixed points. 

6. Charge quantisation on rotational domains. 

We now have expressions for the modes and heat-kernel on S^/F in the ro- 
tational case and have tacitly assumed that all integer values of 2q are allowed. 
However 2q need only form a subset of the integers. A geometric, Dirac type argu- 
ment will firstly be used to calculate this subset. 

Choose a fundamental domain JF which does not contain the string. Let L G J-' 
be a piecewise continuous loop in JF and construct the parallel propagator, 

J[L] = e""/^^ = e-'^", (38) 
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where A is the potential one-form, O is the area within L. This foUows from the 
exphcit form of the monopole potential. As L shrinks to a point in JF, T[L] clearly 
tends to unity as the area vanishes. Now L also describes a loop with (oriented) 
area — (|JF| — O) in S^/F, since S^/F is closed for rotational F (except for the fixed 
points which we can ignore for these purposes as a set of measure zero) . Thus if we 
let L expand to fill the boundary dJ-' C S^, the propagator will also take the value 
1. This requires that q\J-'\ be a multiple of 2n which gives the possible values of the 
magnetic charge as 

Q = |F|-, riGZ, (39) 

the expected value since the magnetic charge per closed domain, q = q/\r\, equals 
the Dirac value, n/2. 

If the string passes through JF there is an extra phase —Arrq in (38) which makes 
no difference to the final result. 

For the cyclic group, Cf., according to (33), and (36), 

a(7*P) = e^-^P^/" = (-1)"P (40) 

so all the 0(7) are ±1. 

A further restriction occurs if Cj^ is a subgroup of a point group F, for then 



a(7-) = e^- 



r\n/k 
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but, by observation, for a noncyclic point group, all |F|//c are even and the 0(7) 
are unity. Furthermore the monopole charge q is integral and there are no spinor 
modes. 

The quantisation (39) or (36), will now be derived in another way. The con- 
sistency condition (32) must hold for both fixed points of 7. Geometry shows, if 
orientation effects are taken correctly into account, that A^(— r^) = — A^(r^). This 
also follows from the mode transformation (31) and the behaviour under parity 
(14), which involves q -^ — g. 

The two definitions of 0(7) require 

k 
Q = -^n^ (41) 
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for some integer n^, where k^ is the order of the cyclic subgroup generated by 7. 
Since A^(r^) = 27r/n^, all the 0(7) are ±1. 

If F is not cyclic, the minimum requirement of (41) is q = {K/2)p where p G Z 
and K is the LCM of all the n . Looking at the character tables of the double 
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point groups it is easily checked that it is not possible to reconcile the tabulated 
±1 representations 0(7*) with (41) unless all the a(7*) are unity, implying that p is 
even or q = {2K/2)p' for p' G Z. But 2K is nothing but the group order F leading 
to (39) for all n. QED. 

7. Integrated kernels and zeta functions 

Special interest attaches itself to the integrated kernel, 

KAr)= Ksvr(r,r;r)rfr = Va(£^) / {T iKs2){r,r;T)dr 

^rn5Z^(^)/ {Tg-iKs2){r,r;r)dr. 
I I ggr -^s 

It is convenient to write everything in terms of the covering S^ quantities because 
we can use the mode form (21) in (42) to obtain 

KAr) = Y. Ml)e-'^^ (43) 



with 

r 



Ml) = ^^J2''i9)x^'^i9-'), (44) 

g 

where we have used the orthogonality of the S^ monopole modes and x (q) is the 
usual SO* (3) ~ SU(2) character of g. In particular 

/ TY [(T^Y^,)(r)Yj,(r)] dr = Xiig)- (45) 

Algebraically we can see that the degeneracy, d^il), is the number of times the 
irreducible representation 'a' occurs in the /-representation. Since degeneracies are 
real, a{g) in (44) and in the following equations, can be replaced by its real part, 
Kea{g), although in the present case the minimal a{g) = ±1. 

(44) shows that the summand in dj-{l), is a class function so that we can make a 

convenient geometrical decomposition of the traced kernel, as in our earlier work [1] . 

The preimage sum, i.e. the sum over F, is firstly replaced by a sum over conjugacy 

classes, {g}, 

2/ _|_ 1 1 

' ' ' ' {g} 
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with g being a typical element in {g}. 

We now recall that the elements of a class correspond to rotations through 
one fixed angle about a set of conjugate axes. For a given set of such axes, one 
corresponding class can be considered to be the primitive class, all others associated 
with these axes then being generated by this one. Thus the sum over all classes 
can be rewritten as a sum over primitive classes and powers of these. Let k be the 
generic order of the rotation associated with the generic primitive class {g} so that 
g^ = id. Then |{^}| is just the number, n^,, of conjugate /c-fold axes and we can 
write 

dAl) = -^ + ^E^^E«'(^)^^'^(^') (46) 

which is the same as when there is no magnetic field, apart from the restriction 
I > \q\ and one sees again, c/ , that the cyclic groups, C^,, form the basic building 
blocks. This can be made explicit as follows. In the case that F is just Cj^, one has 
from (46) 

^^(^) = ^E«'(^)^^'^(^) = ^ + ^E«'(^)^^'^(^) (47) 

p=0 p=l 

so that, substituting back, 

dAl) = ^(1 - E^^) + ^ E^^^^^(0. (48) 

9 9 

where di{l) = 2/ + 1 is the full sphere degeneracy. 

Having now the degeneracies and the eigenvalues, we can turn to the explicit 
construction of the integrated heat-kernel (43) and its Mellin transform, the (- 
function, 

Cr(^) = E^- (49) 

I ^i 



The eigenvalues are given by (2) and we face the old problem of computing 
spherical spectral quantities. We will approach this by firstly looking at a 'linearised' 
system, one whose eigenvalues are 1 + 1/2. The corresponding heat-kernel, denoted 

'^ 1/2 

K-p{t), can be considered as that for the pseudo-operator (-^32 -|- 1/4 -|- g^) and 
will allow us to find the C^function for the eigenvalues (/ -|- 1/2)^ quickly and that 
for the A; of (2), more elaborately. The expressions are also of some statistical 
mechanical interest if r is interpreted as an inverse temperature. 

14 



From (43), we have the hnearised, or square root, kernel 

oo 

KAr) = Y. rfr(/)e-('+i/2)- (50) 

with degeneracies given by (44). 

In accordance with (48), it is sufficient to consider the cychc group case, Kj^{t) 
and Cfc(s)- 

Since the cyclic axis is immaterial, we choose it to lie along the ^-axis and could 
use the expressions in the previous section since df.{l) is TrP, where P is given by 
(35) with r = q. Therefore 

^^ 1 °" ^ 1 p2'Ki{m-\-q) 

^ki'^) ^ ^Z^ Z^ ^ ^ l_e27Ti{m + q)/k- (^1) 

l = \q\ m= — l 

The summations can be relabelled using (5) and performed, but it is perhaps more 
elegant to substitute (44) into (43) and do the /-summation first, as in [1]. 

Remembering the charge quantisation condition, (40), the degeneracies are 

..a..)4gco.(.^) -'""^';y . ,^n./, (52) 

k ^ sm{Tvp/k) 

where n is even or odd. If n is even, the degeneracies are identical to the case when 
q = 0, df,{l,q) =4(/,0). 

We now introduce the generating function 

oo 

h,{a,q)=Y,dk{l.Q)^' (53) 

Mil 

closely connected with the traced heat-kernel, (43), if a = e~^. For n even the only 

effect of the monopole is to make the series start at / = \q\. However it is better to 

continue with the summations. We have, for both integral and half odd-integral q, 

u r \ ^ "^ ST- I . sin(2/ + l)7rp/A;) ; 
^^(-'^) = ^ E i:cos(.rnp) a 

Z=|q|P=0 ^ ^1 > 

_ 1 v^ / N v-^ sin(2/ + l)7rp//c) ; 

k ^-^ ^-^ sin(7r»//c) 

P=0 Z=|g| ^ ^1 > (54) 

(l + a + 2g(l-a)) 1 ^^ ^ sin(2/ + l)np/k) 



^ ' P=l l=\q\ ^ ^' ' 



1 1 + (7^= 2q 

+ 



\ — a \ — a^ ' k{l — a 
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and so the corresponding heat-kernel is related in the simple way 

KI{t) = e-^- (^K,{r) + ^K^{r)^ (55) 

to the monopole-less {q = 0) expression, [1], 

~ ^ coth(fcr/2) 
^^ ' 2sinh(r/2)' 

The decomposition into conjugacy classes, (48), shows that this relation will 
follow through for all groups F, 

i^^(r)=e-^-(KrM+2#..(r)^, (56) 

where g is the charge per domain and Ky{t^ has been determined in [1], 

k frl = cosh(rfor/2) 

^^ ' 2sinh(c/ir/2)sinh(rf2r/2)' ^ ' 

Here d^^ d^ and d^ are integer invariants associated with the reflection group having 
r as its rotation subgroup. 

The zeta function corresponding to the linear heat kernel is easily determined 
as the Mellin transform of (57). In fact we shall find it more useful to consider the 
slightly more general zeta function defined by. 



1 ("^ ~ 



(58) 
= C,2{s,a^q\d^,d<2) + C2('S, a + 9 + c^oMi. '^2) + ^ChI^, a + q), 

where C,ij is the Hurwitz zeta function, and C2 is the two-dimensional Barnes zeta 
function defined for s > 2 by, [19], 

00 ^ 
(2(5,01^1,^2) = J] ^ -— -^. (59) 

Cj7 is actually a one-dimensional Barnes C^function. 

To conclude this section we shall present some properties of the Barnes zeta 
function needed later. 

The function C2(s, a|(ii^, ^2) has simple poles at s = 1, 2 whose residues can be 
written in terms of generalised Bernoulli polynomials 

Res 5^^C2(«, a^i, ^2) = -j-^S^-r(aMi, ^2) (60) 

16 



for r = 1,2. Here we have used the more standard notation as in Erdelyi [20]. 
The values of the Barnes C^function at negative integers are also given in terms of 
generalised Bernoulli polynomials. For n E Z we have, 

Cof— n, a|(ii , dn) = — —B^2„(a\d-,,dn). (61) 

The explicit generalised Bernoulli polynomials required in this paper are, 

So (aMi, (^2) ='^^ ^ (d-i + d2)a + -7(^1 + 'idid2 + d^) 

6 
01 1 

S3 (a|(i]^, ^2) =a"^ — -^{di + 6^2)0^ + -(^i + 'idid2 + ^2)'^ ~ t('^i'^2 + '^i'^2)- 

Finally we present a useful Bernoulli identity which will be used to simplify some 
expressions later on, 

5i'n^i + ^2-«Mi,^2) = (-l)"5i2)(a|rfi,rf2). (63) 



8. Extension to reflection groups 

Before showing how to deal with the eigenvalues (2), we extend the analysis 
to orbifolds S^/F' where F' is a finite reflection group - the complete symmetry 
group of a regular solid as outlined in section 1. The domain of interest is a Mobius 
triangle on S^. 

As shown in section 3, under reflection, the magnetic charge, g, changes sign 
and the projection has to take this into account. The rotational projection is given 
by (22) and (27) and all that is necessary is to combine this with the group decom- 
position (1). We start by writing down the projected modes 

W^ii^) = lS')(r) + a(a)T^y_^S(r), (64) 

where o" is a reflection, say in one of the symmetry planes, so that a{a) = ±1. We 
can choose either sign. 

Note the extended periodicity condition 

T^,M^« (r) = a(7)a(a)VF«^(r), 7' = 7a. 
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which shows that the wavefunction on the domain 7'jF has monopole charge —q if 
that on JF has charge q. 

The group decomposition of V gives 

= U7(l + ^)^, 

and we see that the only possible theory using images requires that adjacent domains 
on the sphere have opposite numerical monopole charge. The charge is q on 7JF 
and —q on 7'JF for all 7 G F, 7' G F'^. 

For the theory to be consistent, we must at least show that the fundamental 
domains with charge q and —q define equivalent physical theories. There are two 
points. The first is that the sign of the charge is arbitrary, being essentially a matter 
of definition for the observer. The second point is that any physically significant 
quantities will depend only on F^^F^^ ~ q^ . Our theory therefore has the possibility 
of being consistent and we now derive the values of q for which it is consistent. 

Across the boundaries of the fundamental domains we have B — i> — B and 
A -^ —A. For consistency we should define these vector quantities to vanish on 
the reflecting boundaries. Consider now the parallel propagator, (38), where the 
loop L ^ T . Since we have defined A to vanish on dT we have X\dT\ = 1. Thus 
for an arbitrary loop approaching the boundary we require trivial parallel transport 
just as in the pure rotational case. This gives the possible values of q as 

n' 

g= |F'|— = |F|n', n' G Z, (65) 

so q is integral and there are no spinor modes. 

The existence of fixed points imposes certain conditions. The situation is more 
restricting than in the pure rotational case because the fixed points form a contin- 
uous set, the boundary of our domain, JF. Let a be a refiection (a^ = 1). The fixed 
point set, V, is the intersection of the refiecting plane with the S^ i.e. a great circle. 
Extending the periodicity condition (18) to F' we see that ifj would have to satisfy 

e-''i^-'^''^+"i''ij{r)=a{a)ilj{r), Vr G P. (66) 

Since ficr(r) is not constant on V, this equation cannot be satisfied on all ofV unless 
one of g = 0, V'lp = 0, n ± P or n G P is true. Applying this argument to the two, 
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or three, reflection generators removes the last two possibilities and we are left with 
either g = or f/'lp = (or both). 

The fact that the magnetic charge, q, has to be of opposite sign on adjacent 
domains under the action of V for the image method to work indicates, crudely, 
that q is zero on the boundary and suggests that our construction satisfies the 
restrictions following from (66). 

Developing this nonrigorous analysis, since the magnetic field B vanishes on 
the refiecting boundaries we must insist that the monopole modes take the value 
Wq^ on the boundaries. Thus the consistency equation (66) is satisfied. Let 7 be 
an arbitary rotation in F. The rotational consistency (32) is satisfied for all 7 since 
q is an integer multiple of |r| and, as already stated, 0(7) = 1,V7. It would thus 
appear that all the values of q in (65) produce a consistent theory. 

For clarity we restate that the mode labelled by qlm takes the values 

W"i2 on FJ^ 



qm 



W'_^>^ on F;^ (67) 

wl^l on a(F^)~-a(F;^). 

The modes (64) can be used to define a heat-kernel analogous to the rotational 
case. We skip straight to the linear heat-kernel which may be written 

» 00 

K^,{r)= / Ve-(^+i/2)-TV[W«(r)W«t(r)]rfr. 

Jsyr' ,^^ 

Here Wg is the vector of solutions Wqm. We can extend this integral to all 
S^ just as we did in equation (42), in the rotational case (this is not entirely trivial 
but it is possible using the invariance of the theory under q -^ —q) , 

^r'ir) = ^J f;e-('+V^)-TY[w(0(r)W«t(r)]dr. 

L=q 

The next step is to the use the explicit rotational modes (27) or (29) (with a{g) = 1), 
equations (13) and (17), and the invariance of the heat kernel under charge reversal. 
The result is, 

K'Ar) = '-K^Ar) + ^ E / Ee'^^+^/^^^Tr [(T^,Y« )(r)Y(Ot(r)]rfr. (68) 
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This equation actually represents two heat kernels for the cases a{a) = ±1. 
For a{a) = —1 we shall use the term 'Dirichlet' and write K'j^ij). For a{(j) = +1 
we use the term 'Neumann' and write K'}^{r). These names are used in analogy 
to the case with no monopole field, g = 0. In this case W^j^ is, by construction, 
a solution on the whole sphere (see (67)) and must satisfy T^/VFq^ = a(a)VFom 
everywhere. On the reflecting boundaries, a{a) = —1 requires WqI^ to vanish i.e. 
Dirichlet boundary conditions, and a{a) = +1 requires the normal derivative of 
Wq^ to vanish, i.e. Neumann boundary conditions. 

We now turn to some speciflc calculations of the heat kernels. Using (45) we 
can write the second term in (68) as. 



a{a) 



/ E E X^(7>-(^+^/^^^ (69) 



'S~ l=q ^/gr^ 



Our flrst calculation is for the reflection group F' with rotational subgroup Cj^. In 
this case F'^^ consists of k reflection planes with a common invariant axis, and with 
angle 2tx /k between adjacent planes. If we take one plane to be the z — x plane then 
we can write F'^^ = {Iig^'^P\p = 0, 1, 2, . . . , /c — 1} where 11 is parity, g^ is a rotation 
by angle tt about the y axis, and 7 is a rotation by angle 2ti /k about the ^-axis. 
Using the explicit result 

V^'^^^,{g^) = {-iy-^d^^, (70) 

and equations (15), (34) we flnd ^^(7') = 1 for all 7' G F'^^ and all I. Thus (69) is 
trivial to calculate in this case, and from (68), (55) we find 

_ -fcr/2 -qr 

K'^ (t) = e~'^'^ - [- — - ('711 

^^ ^ 4sinh(r/2)sinh(A;r/2) 2A;sinh(r/2) ^ ^ 

_ A;t/2 p-qr 

R'5 Crl = p"*?^ - I ^ frol 

^^^ 4sinh(r/2)sinh(A;r/2) 2A;sinh(r/2)' ^ ' 

The first term in each of these expressions is simply exp(— gr) times the monopole- 
less heat kernel [1]. Notice that the extra monopole contribution is the same in 
both heat kernels. 

We can also calculate the heat kernels for the group F' with dihedral rotational 
subgroup Df.. We take F'^ as crDf, where cr is a refiection in the x — y plane, and Dp, 
is the dihedral group with the k-iold cyclic group about the z axis, and a rotation 
by n about the y axis. Using (34) and (70), we find XiiY) = 1 ^or elements 7' G F'j^ 
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involving the rotation by angle tx about y. These elements are |r| in number and 
contribute to (69) the expression 

^^I!^ (73) 

8 sinh(r/2)' ^ ' 

The remaining elements can be written Hg^^^ where g^ is a rotation by tt about the 
z axis, and H, 7 are as defined above. Using (15) and (34) we find the transformation 
matrix (— 1)'+"^5^, exp(—i2'n'mp/ k) . This transformation is non-trivial and we have 
to construct a sum similar to (51). 

For technical variety we use the relations (5) so that the sum can be written, 



CXD 00 



a{a) 

-^"^ ' z^ z^ ^-^^ 

u=0 t= — 00 

= a(a)e-^ 



u=0 t= — 00 

cosh(A;r/2) 



(74) 



8cosh(r/2)sinh(/cr/2) 
Adding together (73) and (74), and using (55) gives the results from (68), 

-(A:+l)r/2 -qr 

K'^ (t) = P~'l'^ - I ^ ^ (7K) 

^^ ' 4sinh(2r/2)sinh(A;r/2) 4)fcsinh(r/2) ^ ' 

_ „(A;+1)t/2 -qr 

Z-q ( ) ^ p-qr "^ L _y "^ (n(^\ 

^^ ^ 4sinh(2r/2)sinh(A;r/2) 4fcsinh(r/2)' ^ ' 

Again the first terms are simply exp(— q'r) times the q = case. 

For all the heat kernels derived above the extra monopole contribution is simply 

given by 

n f,-qT 

(77) 



q e-1^ 



\T'\ sinh(r/2) 

Since we may construct the heat kernel for an arbitrary reflection group from the 
heat kernels calculated above, [2], the simplicity of (77) leads us to conclude that 
in the general case, 

~ p-doT/2 -qr 

K^ Crl = e~'i'^ - V- n' - ('781 

^^ ' 4sinh(rfir/2)sinh(d2r/2) ^ sinh(r/2) ^ ' 

_ pdoT/2 -qr 

^^ ^ 4sinh(rfir/2)sinh(rf2r/2) ^ sinh(r/2)' ^ ^ 

These two equations are the culmination of this section. 
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The first terms in these heat kernels are just exp(— q'r) times the q = expres- 
sions [1], and we have defined the monopole charge per refiection domain 



ir' 



5 m,| 



Equation (68) predicts that the linear Dirichlet and Neumann heat kernels should 
add up to give the rotational linear heat kernel (56), and this is seen to be true (in 
performing this sum we must use the relation q = 2q' for the same value of q). 

The linear zeta functions for K"^ and K^^ are calculated using equation (58) 
with K^{t) suitably replaced. In terms of the Barnes zeta function we find, 

(i^(s, a) = (2(5, a + q + d^ld^, d^) + q'Cnis, a + q) . (80) 

Q(s, a) = (2(5, a + q\d^, d^) + q'Cni^^ a + q) (81) 

It should be remembered that although we may add these two zeta functions to 
produce the rotational zeta function, the physical theories are completely different 
due to the changing sign of the magnetic field. Thus the procedure can only be 
regarded as a formal trick. 

9. General zeta function and its derivative 

Equations (58), (80) and (81) give zeta functions for the eigenvalues (/ + a)^ if 
s is replaced with 2s. In general we need the zeta functions for the more general 
eigenvalues A^ = (/ + a)^ — a^, with suitable constants a and a. For example we 
may add curvature coupling and mass terms to the Hamiltonian Hg2. In this case 
the eigenvalue equation is, 

(iys2+ei? + m2)y« = A,y« 

and we find a = 1/2, a^ = g^ + (1/4 — 2^) — m^ (R = 2 on the unit two-sphere). It is 
assumed that a^ is positive. To analyse these general zeta functions we use similar 
methods to those found in [21]. For brevity we shall just write ({s) to represent a 
general zeta function. For s > 1, we have the explicit mode sum 

as) = E Tn ?o^ 21. - (82) 
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The function d{l) is the degeneracy of the eigenvalue A; for some hnear zeta function 
(i.e. rotational, Dirichlet or Neumann). If we assume that |a| < \q\ + a we can 
perform a binomial expansion on the summand which leads to a continuation of 
C(s) given by 

r=0 ■ ^ ' 

In the the above equation ('^{s, a) is intended to be any one of the rotational, Dirich- 
let, or Neumann zeta functions defined by equations (58), (80) and (81) respectively. 
In order to tie in with [21] we generalise to the case where C^is) represents 
an arbitrary zeta function on a d-dimensional space with simple poles (only) at 
s = 1, 2, . . . , (i. This is the situation encountered in [21] where C'^{s) would just be 
a (i-dimensional Barnes function C^(s,a|d) (of course the label q is defunct in the 
general case). Near the poles we define, 

for r = 1, 2, . . . , (i. For our three cases (with d = 2) the residues A^^ for r = 1, 2 can 
be calculated from the specific forms of the zeta functions and (60). 

The important fact is that the series (83) reduces to a finite sum when s is a 
negative integer. Thus we concentrate on these values of s and find for n G Z"*", 

(-a^)-"C(-n) =5:(-a^)-'(")c'(-2r,a) + l5:a-!li!-lil:iV,,. (84) 

r=0 r=l ^ ^' 

The number u in the above is defined as [d/2] where d is the dimension of the space 
under consideration. The derivative of the zeta function at s = —n can also be 
calculated from (83), 

n 

(-a')-X{-n) =2 5^ (^^)(_«2)-^.'(_2r,a)- 

n 

'^° I (85) 

E«"^^^^^ 1^2.+ 2^2. Wr) -^(n + l))| + 

i!(r — 1) 
(r + n)l 



r=l 

oo 



E "^'fi^c'(2.<.) 



r=u-\-l 
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where i(^{z) = r'{z)/r{z) is the logarithmic derivative of the gamma function. 

The problem now rests on the evaluation of the infinite sum on the last line of 
this expression. We show that this sum can be written in finite terms. It is expected 
that the sum will be finite since we expect ('(—n) to be finite, and all other terms 
on the right hand side are finite. A word of caution is that in singular situations, 
there is the possibility that logarithmic terms, logr, may appear in the asymptotic 
expansion of the heat-kernel. If this were so, more care would have to be taken over 
the evaluation of the determinants. However no such terms occur here. 

Using the integral representation (58) extended to the arbitrary zeta function 
C^is, a), we may write the last line in (85) as, 

2^' E «' ^ ^ ;„o „ / T^'-'e'l^-''^K\T)dT, (86) 

n=u+l ^ ^ ^ ' •^^' 

where K'^{t) is the linear heat kernel associated with C'^{s, a). Since (86) is assumed 
finite we may take the sum inside the integral. Thus our problem can be reduced 
to evaluating the sequence of sums. 



Tnir) = nlJ2 



j.\r^2r 



-^ (2r)!(r + n)!' 

Using the simple result y^(2r)\ = 2^^r!r(r + l/2) and changing the summation 
variable to r' = r + n gives the result 



-r^l 



1 

T^{r) = n\V^{-r)'/'-I_^_y,{r) - J^ (^)(2r)!(- 

r=0 

where I^{x) is the modified Bessel function. 

To find a closed form for T^{t) we employ a useful integral representation given 
in reference [22]. For z/ > 0, 

r(i + z.)/_(x) = 



2 1 



e-"*(l - t^y-^/^dt + sin(7rz/) / e-"*(t2 - If-^'^dt 



Setting z/ = n + 1/2 in this expression we may expand the integrand factors (1 — t^)" 
using the binomial theorem leaving simple exponential integrals. After a little work 
we find, 

Tn{r) = J2 O i-^y {{-r')-''{2rn2n - r)!^ (e^ + (-l)'^e--) - {2r)lr-'A . 

r=0 ^ ^ 

(87) 
24 



Having found a suitable expression for T^ we can now go back to (86) and write 
it in the new form, 

r jarja.) - ±(--f \r"ln2Ziy] r-'^"'-"I<^(r)ir. (88) 

We have introduced into the integral a regulator r* (the expression that we want is 
given for the value s = 0) . The continuation variable s has been introduced so that 
we may evaluate the integrals of the individual terms in the sum definition (87) of 
T„, and of the sum subtracted from it, before we perform the sums. We assume 
that s is large enough so that all the individual integrals are well defined. In fact 
this requires s > 2n. Performing the integrations leaves, 

n 

(-«')"" Y. C") (2«)'(2^ - rV.r{s + r - 2n) X 

{(''{s + r-2n,a + a) + {-lyCHs + r - 2n, a - a)} - 

2i:(-°")-l;)(2'-)!n»-2r)C'(»-2r.a)-X:a'- "-;.'--y';|.'_+y C'(»+2n^ 

(89) 
As s -^ all of these terms diverge, although taking all terms together we must get 
a finite result i.e. all poles must cancel as s ^ 0. This cancellation of the poles 
leads to the equation. 



■^7 

n 



(-«')"" E (^) (2«)' {C(r -2n,a-a) + {-iyC{r -2n,a + a)}- 



.r 

r=0 



2E(-"^)-i:)C'(-2''.a)-5:""^fe#'V,, = 0. (90) 

Comparing equations (90) and (84) we see that the pole cancellation is precisely 
the statement, 

C(-^) = 2 E (^) (2«)'^ {C(^ - 2n, a - a) + {-iyC{r -2n,a + a)}. (91) 

r=0 

This expression contains equally terms with arguments a + a and a — a. We shall 
say that C(^^) is 'symmetric'. It is a generalisation to general n of the symmetric 
expression for n = found in [21]. The methods used in this reference do not 
produce a suitable pole cancellation to give a symmetric result for C(— n). 
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The finite remainder part of (89) as s ^ is given by tfie expression, 

n 

(-"')"" Y. C") (2«)' {C'{r -2n,a-a) + {-lYC'{r -2n,a + a)} 



r=0 



2^n!(r-l)! 



n u 

r=0 r=l ^ ' 



Ror + 



n 

(-a^)-" Y^ (^ V2a)X2?i-r + 1) {C{r - 2?i, a - a) + (-l)'^C'^(r - 2n, a + a)} 



r=0 



r=l 
n-1 



Inserting tfiis expression into (85), and adding zero in tfie form of ?/;(2n + 1) times 
(90) gives 

n 

C'(-^) = E (^) (2«)'^ {C'(^ _ 2n, a - a) + (-l)^C'(r - 2n, a + a)} - 

n 

E (^) (2tt) V, {C(r - 2n, a - a) + (-1)"C'(^ - 2n, a + a)} - 

r=l 

E(-a^)''"'^ (^) (2V'(2r + 1) - V(r + 1) + ^(n + 1) - 2^(2?i + l))C'^(-2r, a)- 

(-1)- E a^(-+") ^7 ~ ,, (2V'(2r) - V(r) + ^(n + 1) - 2^(2n + l))N^,. (92) 

The quantities a^ in the above equation are defined as the sums, 

r-l ^ 

'^'^"■^2^r^' 

A;=0 

Equation (92) is not symmetric in the sense that, unhke (91), it does not depend 
only on the quantities a ± a. If we assume that the sum over the residues is a true 
feature of C,'{—n), as it is for C'(0) i^i [21], then we are stiU left with a sum over 
(■•?(— 2r, a). We will now re-write this sum in a more natural, i.e. symmetric, form. 
To this end we introduce the intermediate zeta function, C,{s), on the space 
JA = M^" X S^/r, which is given by 

lU) = il^_!^^(s - n) . (93) 

^^ ' (47r)"r(s)^^ ' ^ ' 
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Combining this with (83) gives an expansion for ({s) 



oo 



as) = J2a'^^^^^r{2s + 2r^a) (94) 

n r\T{s) 

where we have defined the new hnear zeta function via 

r(«,a) = 7^%^Cn«-2n,a) (95) 

(47r)"r(is) 

The dimension of A^ is d = 2n + d (remember, for our monopole case d = 2). We 
see from (95) that ( (s) has poles at s = 2, 4, . . . , 2n and s = 2n + 1, 2?i + 2, . . . , (i. 
The residues are given by the formulae, 

^2.= tA ^nr^^~iw ^C^(2r-2n,a), r = l,2,...,n (96) 

{An)^[r — l)\[n — r)i 



r(^r) 

(47r)"r(ir + n 



N,^.^ = ——^f^^—N,, r = l,2,...,rf. (97) 



The purpose of making these new definitions is the functional similarity between 
(83) and (94). This implies that results like (91) and (92) should exist for C(s) in 
terms of ( (s). The important point is that we know that C'(0) can be written as a 
symmetric part and a sum over the residues of C'^{s, a), either from [21], or setting 
n = in (92). Thus we might expect C (0) to consist of a symmetric part and a 
sum over the residues of ( (s, a). Differentiating (93) and setting s = gives, 

C'(0) = -^-^^ iCi-n) + mn + 1) - V(l))C(-n)) • (98) 

Now C(~^) is symmetric in terms of C,^{s, a ± a) and extends easily to a symmetric 
form for C, (s, a) using (95). Thus by our reasoning we expect C,'{~n) to contain a 
sum over the residues A^^. This is exactly what we find, and the final result, in this 
section, is the symmetrical expression, [23] , 



n 

C'(-^) = E (iJ) (2«)'^ {C^'(^ _ 2n, a - a) + {-lYC'ir - 2n, a + a)} - 

r=0 
n 

Y, (^) (2«) V, {C(r - 2n, a - a) + {-iyC{r ~2n,a + a)}- (99) 

u 2r 



2r: 

r 

r=l 



27 



with the definitions u = [d/2] and 

1 1 

p^ =ip{2r - 2n + 1) - -ip{r - n + 1) - (V'(2?i + 1) - -^(n + 1)) 

^^2k + l ^^2k + l 

k=0 fc=0 

The conclusion of these manipulations is that, despite the apparent awkward- 
ness of the binomial expansion, (83), to obtain the required C^function, the quan- 
tities that we want are given in finite terms, (91), (92), (99), and involve only 
relatively standard functions such as generalised Bernoulli polynomials introduced 
via the Barnes ^-function. 

10. Vacuum energy calculations 

Simply as an example of the use of the preceding expressions, we evaluate some 
vacuum (Casimir) energies on S^/F'. 

Let C'^{s,a) represent one of the rotational, Dirichlet or Neumann linear zeta 
functions as in the previous section. Then ('^{s,a) can be extended to the odd- 
dimensional space-time M x S^/F (or M x S^/F') by defining a new zeta function 
({s) given by (with a = 1/2), 

This zeta function corresponds to the rather artificial case of a conformally coupled 
field in three dimensions with mass q^. f The vacuum energy associated with this 
physical situation is defined by the simple formula [24], 

2s 



--i-si^i «" 



(101) 



s=0 

In this equation p is an arbitrary mass scale and fu,^ = 1/r is the mass scale associ- 
ated with the sphere radius r (which has the value r = 1). Equation (101) is in fact 
just half the logarithmic determinant. Inserting equation (100) into (101) gives the 
simpler expression 

£;=^C(-l,l/2), (102) 



jFor a scalar field conformally coupled in A^ dimensions we have 4^ = (A^ — 2)/(N — 1). In fact 
we only require that a^ = or equivalently m^ = q'^ + (1/4 — 2^). 
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which is finite and independent of fu. (we have set r = 1 again). 

The vacuum energy on the space R x R "^ x 5'^/r(') can also be found and 
gives E proportional to C^i—l — 2?i, ^). The calculation for general n is entirely 
equivalent to the n = case, which we now calculate. 

Using the definitions (80), (80) and equation (61), we find from (102) the 
Dirichlet and Neumann vacuum energies. 



-"J 



E{oy = ±^^{dl -dl- 4) + ^{3dl -di-dl + l± 6doq - 2q'). 

The constant (q independent) terms are exactly the same as those calculated in 
[1] for g = 0, as required. Adding together the Dirichlet and Neumann vacuum 
energies gives the rotational vacuum energy 

Er = ^{Ml-dl-dl-2q'), 

where we have used the relation q = 2q' for fixed q. This vacuum energy necessarily 
vanishes for g = 0, as proved in [1]. We now list the vacuum energies Er d\ for all 
possible refiection groups F', 

29 a' 
O*: ±— + 1^(383 ±108q-V) 

89 a' 
y*-- ±^ + 1^(1079 ±90q-V) 

Dn]Cn ■■ T^ + ^(4n2 -l±6nq- 4q') 

D^ineven), D.jDjn odd) : ±^^^±^l|^^^|(4(n + l)(n + 2) - 9± 

192n 48 

6(71+1)^-4^2). 

Those for the corresponding rotational subgroups are obtained by adding the D and 
N values. 

We now go on to calculate the vacuum energy on the even- dimensional space 
r2" X 5'2/r('). Using equations (93), (98) and (101) we find the expression. 



^ = S^ I ^'<-"' ^ 



ln(-^] +^(n + l)-^(l) 



a-n) . (103) 



Here ({s) is the general zeta function defined via equations (82) and (83). 
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An infinite contribution has been (arbitrarily) dropped to arrive at (103), the 
logarithmic term being a relic of this divergence. Since ({0) is not zero for any 
of the three monopole theories, we conclude from (103) that the vacuum energy 
is explicitly dependent on the arbitrary scale fi. However for simplicity we shall 
assume n = fi^ = 1 for the rest of this section. Our concern in this paper is not 
with realistic quantum field theory considerations. 

Equations (91), (92) imply that the calculation for increasing n merely requires 
the evaluation of more and more zeta functions and their derivatives. Thus we shall 
concentrate on the simplest case n = corresponding to S^/r(') itself. This will 
also allow us to compare with the results for g = studied in [21]. 

We consider the case of a massless field with minimal coupling, that is m'^ = 
and ^ = 0. Thus we have a = 1/2 and 2a = \/ ^q^ + 1 (we shall still write a when 
convenient). From the definitions (80), (81), and equations (91), (61), we find for 
Dirichlet and Neumann boundary conditions 

^{"}^°^ = ]^ + ^(^o(^o - 1) + 1 ± 6rfo9 + 18q'). 

These expressions reduce to those found earlier for g = [21]. For the rotational 
case we simply add the Dirichlet and Neumann results to give 

^r(O) = ^ + 7^(^o(^o - 1) + 1 + 18?'). 

D D|i I 

To calculate the vacuum energy we have still to calculate C'(0)- Considering 
equation (92) with n = requires the evaluation of the residue A^2- -^o^ both 
Dirichlet and Neumann zeta functions we find from (60) and (62) the value N^ — 
2/|r'|. Using the derivative of the Hurwitz zeta function in [20] then gives the zeta 
function derivatives 



Cd(0) =C2(0, 2 + ? - « + ^o) + C2(0, - + ? + « + c/o) + 

q' In (r(^ + q - a)r( J + g + a) 1 - q' ln(27r) "" ^ ^^' 



•2 " ' '2 " > \ -^ y , 2|r 



(104) 



C(o) =C2(o, ^ + ?-«) + C2(o, ^ +? + «)+ 

q' In |r(^ + Q - a)r(^ + g + a) I - ^' ln(2 



1 + V 



(105) 



2|r'| 

(The rotational zeta function derivative is just the sum of these two.) The triangle 
inequality |a;| + \y\ > \/ x^ + y^ {x,y E M.) implies that 1/2 + q > a. The equality 
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is only met for q = 0, and in this case we have to remove the singularity in the first 
term C2(05 1/2 + Q — a|(ii, ^2) of ^^e Dirichlet zeta function as in 21 (we shall do this 
later). For g > all terms in (105) and (104) are well defined. 

There is no known analytic form for the derivatives of the Barnes zeta functions 
appearing in (105), (104) and so we have to calculate them numerically. To do this 
we obviously need a continuation of the Barnes zeta function which is open to easy 
numerical computation. In [21] several efficient continuations are presented which 
are valid for positive integer values of di and ^2- However as we shall show in the 
next section, it is useful to have an expression which is valid for all ^1,(^2 G IR • We 
shall now derive such an expression. 

Our starting point is the Plana sum formula which we display here [25], 

J2f{n)=hf{a) + f{b))+ I f{t)dt+ 

"^" (106) 

. r f{a + it) - f{a - it) - f{b + it) + f{b - it) 

' Jo e'"* - 1 

To be valid f{t) must be an analytic function in the region of the complex t plane 
a < Re (t) < 6, and the integrals must exist. Applying (106) twice to the sum 
definition (59) of the Barnes zeta function gives immediately, for s > 2, 



(2(5, a|rfi, d^) =-rfi ^Ch(s, ^) + ^ (g_i) ^i '^His - 1, ^; 



+ 



M e2"*-l 12 ^^^^'^^^^ '-(a-id^t) ') + 



°" _dt r 1 





— -^ ((a + id^t)^-' -{a- id^t)^-') 1 + 

f^ du r dt (, ., .,^,_, 

/ -7^ / -TT-i < (a + id,u — td'yt) + 

Jo e^-^-lJo e2-*-l\^ ^ "^ 

(a — idiU + id2t)~^ — (a + idiU + id2t)~^ — (a — id^u — id2t)~ 

(107) 
It is simple to verify that we are meeting the conditions required for the validity 
of the sum formula. In order to get rid of the single integrals in (107) we use the 
Plana sum definition of the Hurwitz zeta function which is, from (106) 

Us,a) = -a-^ + -+z/ ^;^^-{(a + zt)-^-(a-zt)-^}. (108) 

/ s — 1 Jo ^ "J- 

The integral part of this expression is equivalent to the integrals appearing in (107). 
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To simplify the double integral in (107) we first perform a change of variables 
from t, u to diU ± ^2^. Following this, we use the easily proved formula, 

{x + iyY + {x — iyY = 2 cos (s tan""*^ y/x) , 

which is valid for Re {x) > 0. After a little work we find the more convenient form 
for the double integral, 

2 /°°,.,^(!^"-(';-°->/°)). (109) 



did^Jo e2™-l (a2 + «;2)s/2 

All the non-trivial d^, ^2 dependence has been absorbed into the function G{w 
which is independent of s and has the explicit form 



G{w) =(e 



2^w .^] I dy 



(e'5i^-l)(e'52(— ^)-l) 
rfy P dy ^ ^ ' 



/o (e'^i2^-l)(e'52("'+?^) -1) io (e'^i(^+?^) -1) (6-^2^-1) 
where we have defined 5^ = 2n/dj^, i = 1,2. This function is symmetric under the 
interchange of di, ^2 as one would expect. The factor (e'^'^'^ — 1) has been included 
into the definition to ensure that G{w) is finite as w ^ 0. 

All the integrals in the definition of G{w) are divergent at their lower limits, 
and the first integral is also divergent at its upper limit. However one can check 
that the combination is well defined. In fact by expanding the integrands at their 
limits of integration we find that for small e > 0, 



27r«; . ^ ) / dy 



G{w) =(e"™ - 1) 

(.Je (^e"i« — ij(^e"2v- w^ — ij 

+ 0(e). 



[e^-^y - l)(e'^2(™-y) - l) 
dy [^ dy 



{e^iv - l){e^2i^+y) - I) J, (e^A'^+y) -l){e^2y -I] 
Thus G{w) is easy to calculate numerically, with the error being of order e. We may 
also make e the lower limit of the integration over w in (109). Since the integrand 
with respect to w is finite at the lower limit, the error incurred will still be 0(e). 

The full expression for our continuation of the Barnes zeta function, after deal- 
ing with both the single and double integrals in (107), is 

C2(s,a|rfi,d2)=-Ta - V 1 2; , 



4 d^d^is - l){s - 2) 2^1^2(5-1) 

IdrCnis, ^) + ^d^^His, ^) + 

d;dj^){''-'^-^'-'^^^^ 

2 f°° G{w)dw cos [s tan~^ (w / a)) 



(111) 



did2Jo e2^™-l (a^ + w^y/^ 
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Although this formula was derived for s > 2, it is actually a continuation to all 
values of (complex) s except at the points s = 1, 2 where there are simple poles. It 
is easy to check that these poles are correct in that their residues match those given 
in (60). From (111) we can calculate the derivative of the Barnes zeta function at 
s = 0, 

C^(0,a|d„d,) = -iln.-i(l-^-^)ln(2.) + iolna(l + i)- 



"l"2 Jo ^ ~ ^ 

(112) 
Here we have used (108) again to convert derivatives of the Hurwitz zeta function 
Cj:^(s, a) at s = —1 into integrals suitable for numerical evaluation. 

The expression (112) can be used directly in equation (105) and (104). Figures 
1 to 4 show plots of i? = — C'(0)/2 for small values of ^. The lines on the graphs for 
T, O, and Y are for labelling only. The graphs for Cf. and Df. are plotted for all k, 
although at the moment we are only concerned with the integral values /c = 1, 2, . . . 
marked with crosses. On each graph the value q of the charge per rotational domain 
is always twice the reflection value q' so as to give the same value of q. 
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Fig.l Vacuum energies for Cf. and Df. with q' = 0. 
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Fig. 2 Vacuum energies for C/. and Dj^ with q' = 1/2. 
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The graphs with q = q = are exactly the same as those calculated in [21]. 
To plot the graphs in this case we have had to remove the divergence in the first 
Barnes zeta function derivative in (105). This is due to a zero mode appearing 
in the spectrum of the operator Hg2. The divergence is logarithmic and can be 
removed by defining the so called F-modular form p2 defined by [19,21], 

lim C2(07 a\dij ^2) = — In a — lnp2- 

The leading term —In a on the right hand side is the divergent zero mode contri- 
bution which must be removed. The form p2 has been calculated in terms of the 
multiple gamma function by Barnes [19,26]. The F-modular form can be more easily 
calculated by using the simply proved relation 

ci a 

(2(5,01^1,^2) =a"'' + <i^\j:^(s,l + — ) + rf^''CiT(■s,l + -l-)+C2(■s,a + <il + rf2Ml,t^2)• 
Differentiating this with respect to s at s = 0, and taking the limit a ^ leaves 
the expression 

- Inp2 = 2 HdM - ln(27r) + C2(0, ^o + l^i, ^^2)- 

For g > 0, the zero mode that we have just removed for q' = is still almost 
zero. This is the reason why the vacuum energy for the Dirichlet zeta function 
increases more rapidly with q' than in the Neumann case. 

11. Theory for a spherical slice 

By a spherical slice of width /3 we mean the space 5"^ = {(6*, (t))\6 G [0, tt] , ^ G 
[0, (3]} with the points (6*, 0) and (6*, 13) identified (here (6*, (j)) are the spherical polar 
coordinates on the sphere which shall be used throughout this section). One might 
also term this space a 'periodic lune'. 

The starting point for the monopole theory on Sp is the solution to the explicit 
differential equation (4) with arbitrary q. As mentioned in the discussion surround- 
ing the differential equation the solutions are characterised by integers u,v, and 
are given by equations (3), (5). We shall now show that it is possible to define 
consistent monopole theories on Sa corresponding to the rotational, Dirichlet and 
Neumann cases already given. The rotational case is considered first. 

Define khy k = 27r//3 and let 7 denote a rotation by angle (3 about the z axis. 
Due to the identification of points in 5"^ at (/> = 0, /?, we will proceed convention- 
ally and take the wave function to be single valued i.e. to have period 13. Also, 
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for rapidity, the monopole charge quantisation wiU be obtained by requiring the 
connecting function exp{i2q(j)) to have period /3. This yields 

k 
q = —n , n E Z . (113) 

Corresponding to these values of q we find 0(7) = (— 1)". If /c is an integer then 
Sf^ is a fundamental domain for the group Cf^ (generated by 7), and the relation 
\Cj^\ = k gives an equivalence between the quantisation conditions (39) and (113). 
Thus we see that the quantisation condition for arbitrary /c is a generalisation of the 
previous theory for C^, (with integer k). From single valuedness, Lp{9,(3) = (/?(6', 0), 
for the monopole harmonics (5), we see from the explicit (p dependence that, 

m = jk — q , jGZ. (114) 

In this more general setting it would appear that we can always choose n to be odd, 
whereas before this was only possible for k odd also. From (5)) we find the possible 
values of /, 

l^u+-{\j\ + \n-j\) , uez+,jez. 

The eigenvalues A appearing in the mode equation (4) are still given by /(/ + 1) — q^. 
As in section 7 we shall calculate the linear heat kernel for the eigenvalues (/ + 1/2). 
Using the values of q, /, and m above gives the sum form for the linear heat 
kernel analogous to (51), with (5), 

00 00 
Klir) = e-^/2 Yl E e-^("+'=(l^'l+l^'+'^l)/2 . (115) 

■U = J = — CXD 

Evaluating this sum gives exactly the result found before, equation (56), except 
that now k is not an integer, and q need not be an integer multiple of k. A corollary 
of this result is that (56) is also valid for the case when F = Cj^ with k and 2q 
odd integers. The zeta function defined by equation (58) is also correct for di = 1, 
3,2 = k and hence we conclude that the vacuum energies calculated in the previous 
section are valid for arbitrary k. 

For the Dirichlet and Neumann theories defined in section 8 we consider the slice 
3^12 and two refiection planes Pq and P^ which leave (/> = and (f) = /3/2 invariant 
respectively. Following the same arguments as in section 8 gives the values of the 
monopole charge 

q = kn', (116) 

39 



which is equivalent to equation (65). Taking 7' in (64) to be a reflection in Pq and 
using (17) gives the modes WqrL as the combinations, 

Ti-W (^, <p) = y«(^, <P) + a(a)y«^(^, </>) 

Here Yqm are those deflned in (3) with the string along the negative z axis, and 
a{a) = ±1 as before (the reflection in Pq is expressed simply as ^ ^ —</>)• 
The (p dependent part of the modes WqJi is given by 



^^ ^ ■' {sm{m(j)) , a((T) = -l 

Reflection in Pq is equivalent to (^ ^ —(p and we see explicitly that this transforms 
WqmiO.O) into a{a)W_q^{6,0) as required. Reflection in P^ is equivalent to ^ ^ 
(3 — (j) and results in the condition on to, 

m = kj , j eZ+ . (117) 

For the Dirichlet case a(o") = — 1 the m = mode is in fact zero and has to be 
removed. Comparing (117) with (114) (with n = 2n' even) we see that the only 
difference between the rotation and reflection cases is that in the reflection case the 
values of m are restricted to positive integers. 

The linear heat kernel for the Neumann case a{a) = +1 is given by 

00 00 
K^(r) = 6-^/25^5^ e-^^^+'^^l^'+^'l + l^-^'D/^. (118) 

u=0 j=0 

Explicit evaluation of the sum gives exactly the heat kernel as before, equation 
(55). The Dirichlet case involves subtracting the j = term from (71) and yields 
the previous expression (75). Thus the zeta functions are given exactly as before 
and we conclude that the vacuum energy calculated for Cj^ is in fact valid for all k. 
To extend the results of Dj^ to arbitrary k we consider the slice of the upper 
hemisphere S'^ = {i6,(j))\e G [0, 7r/2] , (p e [0,/?]} with again k = 27r//3. This is a 
fundamental domain for Df. when k is an integer. The theory then follows as for 
the Cj^ extension above, but we must include in this case a rotation about the x 
axis by angle tt. This rotation can be thought of as a reflection in the plane Pq 
followed by a reflection in the x — y plane, which we call ^2- The theory for 5"^ 
above is adapted to the reflection Pq and thus we see that P2 is the essential extra 
detail here. 
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The reflection P2 is equivalent to the transformation 9 ^ n — 9 which does not 
affect the (p dependence of the modes YgV?i(6', (^). Since the extension to arbitrary 
k is entirely linked with the (p dependence, we conclude that all the heat kernels, 
zeta functions and vacuum energies for Dj^ can be extended to arbitrary k. The 
values of the monopole charge in the reflection case are calculated using the theory 
of section 8 as 

q = 2kn' , n' e Z, 

which is the generalisation of (65) with |r'| = 4k. 

12. Summary and discussion 

We have thoroughly adapted Dirac's monopole theory to the orbifold, S^/F, 
for the cases that F contains only rotations and when F is generated by reflections. 
In the former case we imposed rotational (periodic) boundary conditions on the 
monopole solutions. In the latter we had a choice of boundary conditions defined so 
as to reproduce Dirichlet and Neumann conditions for no monopole charge, q = 0. 
We found that it was the monopole charge q = q/\T\ through S^/F that was Dirac 
quantised with 2q G Z. 

After all the formalities of the theory had been tidied we explicitly calculated 
the vacuum energies on the orbifolds S^/F and M x S^/F. Formal expressions are 
given for the generalisation to the spaces M^"^ x S^/F and M x R^"^ x S^/F. Finally 
we provided an extension of the monopole theory to arbitrary slices of the sphere 
and hemisphere. In this case the flux through the spherical region is still quantised, 
although now the overall monopole charge q is not, in general, an integer or half 
odd-integer. 

We feel that the scalar theory has been developed essentially to its analytical 
limit on the factored sphere. The next step would be the extension to M /F for F a 
reflection or rotation group. This requires modes of the full Hamiltonian which are 
given by 



Ym,<P)Ukr)^/k/^, u = y/{l + l/2y 



q,2 



with eigenvalues /c^. Since the radial dependence does not involve to, the underlying 
facts of the theory (modes on factored space, charge quantisation etc. ) are the 
same as in the spherical case. However the heat-kernel calculation, and hence the 
(■-function, is completely different. Due to the complicated index, a/(/ + 1/2)^ — Q^, 
closed forms do not seem possible and asymptotic methods are needed. One could 
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always arbitrarily add a term (^ jr to the (total) Hamiltonian and then a closed form 
would exist. This fact suggests that there is some significance to this modification. 

The spinor theory on the factored sphere S^/F has been considered by Chang 
[27]. He found a consistent theory only for T = C^ with k odd. For g 7^ we claim 
that the same restriction still holds. This follows from the lack of half odd-integral 
solutions to the scalar monopole problem for V ^ C^^. 

A possible extension of the scalar calculation would be to consider if (high 
temperature) Bose-Einstein condensation occurs. The general theory has been laid 
down by Toms [28] See also Kirsten and Toms [29]. Basically, all that is required 
is to ensure that the ^-function for the theory, and its derivative, are finite at zero 
as the chemical potential approaches a critical value. On the two-sphere we can 
use the calculations of C(^^) ^ind C,\—n given in section 9 to study the theory on 
m}"^ X 32. In fact we could also discuss R^'" x S'^. 
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